Layered Kondo lattice model for quantum critical /3-YbAlB 4 
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We analyze the magnetic and electronic properties of the quantum critical heavy fermion superconductor 
/3-YbAlB4, calculating the Fermi surface and the angular dependence of the extremal orbits relevant to the de 
Haas-van Alphen measurements. Using a combination of the realistic materials modeling and single-ion crystal 
field analysis, we are led to propose a layered Kondo lattice model for this system, in which two dimensional 
boron layers are Kondo coupled via interlay er Yb moments in a J z = ±5/2 state. This model fits the measured 
single ion magnetic susceptibility and predicts a substantial change in the electronic anisotropy as the system is 
pressure-tuned through the quantum critical point. 
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Heavy electron materials have played a long-standing role 
in strongly correlated electron physics, providing key insights 
into emergent quantum mechanical behavior that many hope 
can be scaled up in energy and temperature in actinide and 
transition metal compounds. One of the remarkable aspects 
of heavy fermion materials is that they can be tuned, by ap- 
plying pressure or magnetic field at low temperatures, through 



a quantum critical point 
liquid behavior fl S, i,^ 
wards superconductivity 



,0], where they exhibit non-Fermi 
and a marked predisposition to- 
,0]. The recent discovery of a lay- 
ered Yb based heavy fermion material /3-YbAlB4 [7] which 
is both stoichiometrically quantum critical and superconduct- 
ing at T c = 80 mK has attracted a great deal of interest. 
There are many examples of heavy electron superconductiv- 
ity, often in close proximity to a quantum critical point [Eh- 
However, although heavy electron behavior and quantum crit- 
ical behavior have been observed in Ce, U and Yb materi- 
als Q, H @, EH HH [H [H, /3-YbAlB 4 is the first ytterbium 
based heavy fermion superconductor. This material raises 
many fascinating questions. Why, for example, is this system 
intrinsically quantum critical, with a specific heat coefficient 
that is finite in a magnetic field, but which diverges as the field 
is removed? 

In this paper, we present an analysis of the magnetic and 
electronic properties of /3-YbAlB4. We use a combination of 
electronic structure calculations and crystal field analysis to 
develop a simple model for the low energy physics. Our re- 
sults are consistent with the almost localized nature of the Yb 
f-electrons. In particular, we find that the electron effective 
mass, as measured by the low temperature specific heat coef- 
ficient 7 = C v /T, exceeds the band- structure value by more 
than a factor of 30. Moreover, in an LDA+U band- structure 
calculation, the system is found to develop an antiferromag- 
netically ordered ground-state, underlining the close vicinity 
to magnetic instability. 

We have calculated the Fermi surface and the angular de- 
pendence of the extremal orbits relevant to the de Haas- 
van Alphen measurements. Despite the layered crystal struc- 
ture, hybridization with the /-electrons gives rise to a three- 
dimensional Fermi surface. Combining the band- structure 
calculations and crystal field analysis, we are led to propose 



an anisotropic Kondo lattice description of this material, in 
which Yb magnetic moments at a site of seven-fold symme- 
try hybridize with neighboring boron planes, to give rise to 
the observed 3-dimensional band structure. We show that 
the crystal field configuration with maximum interlay er over- 
lap, | J = 7/2, mj = ±5/2) (Fig. [D, accounts well for the 
anisotropy in the high temperature magnetic susceptibility. 

The compound /3-YbAlB4 crystallizes in an orthorhombic 
structure (group Cmmm) 11411 . with ytterbium and aluminum 
atoms sandwiched between boron layers. In this unusual 
structure the Yb atoms are centered between seven-member 
boron rings, as illustrated in Fig. [TJ while the Al atoms are 
centered between pentagonal rings of boron. A central ques- 
tion raised by the structure, is whether the electronic structure 
will be quasi-two dimensional as in the case of intercalated 
graphite compounds. Our calculations indicate instead that 
the two-dimensional layers of boron are "short-circuited" by 
electrons of Yb and Al. The role of the Kondo effect in this 
setting is particularly intriguing and will be discussed below. 

For the band structure calculations, a full potential linear 




FIG. 1 : (Color online) Local coordination of Yb atom (large spheres, 
blue), suspended between two planes of boron heptagons (small 
green spheres) and surrounded in-plane by a distorted rectangular 
of Al atoms (red spheres). The proposed \mj = ±5/2) ground state 
wavefunction is also shown, see Eq. $4§. Boron conduction electrons 
(Ti±) |0) interact with Yb spins Sj via layered Kondo lattice model 
Eq. (0. 




FIG. 2: (Color online) Projected density of states (DOS), show- 
ing relative contribution of different atoms. The DOS at the Fermi 
level corresponds to the specific heat coefficient 7 = C v /T = 6.7 
mJ/(mol YbxK 2 ) per Yb atom. 

augmented plane wave method (LAPW) was carried out using 
the WIEN2k code LI 5], using a Local Density Approximation 
(LDA) with a generalized gradient correction HI 611 . The strong 
spin-orbit coupling of the Yb atoms was taken into account 
using a scalar-relativistic approach lfl5l [nil. The calculated 
electron density of states (DOS) of /3-YbAlB4 is shown in 
Fig. |2 The Yb /-bands are split by the spin-orbit (SO) in- 
teraction into J = 5/2 and J = 7/2 multiplets separated by 
1 .4 eV, and are heavily hybridized with the lighter conduction 
electrons of B, and (to a lesser extent) Al atoms. The DOS at 
the Fermi level is dominated by the Yb /-bands, with a band 
mass around 5-6 m e in the (afr)-plane and about m e along the 
c-axis. The specific heat coefficient 7th = 6.7mJ/(molxK 2 ), 
extracted from the value of the DOS at Ep contrasts with a 
measured value 7 exp = lim Cm/T > 170 mJ/(mol YbxK 2 ) 

at a field of 0.50 Tesla Experimentally, the zero temper- 
ature specific heat coefficient is field dependent, diverging as 
^-1/4 Q n me b as j s f thi s comparison, 
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where the critical field B c = within current experimental 
resolution. This divergent mass enhancement indicates that 
the Yb /-bands become localized in the zero field limit. 

The LDA calculation predicts that the Yb ions are in a 
mixed valence state, with rif « 13.4. This is clearly incon- 
sistent with the susceptibility measurements, which reveal a 
Curie- Weiss susceptibility at high temperatures, with the full 
moment of the 4/ 13 state: a signature of local moment be- 
havior. We attempted to cure this failure using the LDA+U 
method LI 8] where the Anderson U is treated at the mean-field 
level for the localized Yb 4/-orbitals. These calculations pre- 
dict a magnetic ground state with an in-plane AFM arrange- 
ment of Yb spins and the effective moment of 0.6 (13 per Yb 
ion. While no magnetic ordering is observed in the material, 
we may regard the LDA+U result as a mean-field description 
of strong underlying magnetic correlations. 

Although LDA can not capture the large mass renormal- 
izations of strongly correlated systems, it generally provides 
agood guide to the underlying Fermi surface topology II 19i 
2JJ HI. The calculated LDA Fermi surface /3-YbAlB 4 (see 
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FIG. 3: (Color online) LDA-solution calculated (a) Fermi surfaces 
and (b) angle-dependent de Haas - van Alphen frequencies of ex- 
tremal orbits showing only the lightest bands with frequencies F < 
10 kT (higher frequencies up to 50 kT are also present but are un- 
likely to be observed in dHvA measurement). Two Yb-/ bands con- 
tribute most to the Fermi surface: band 89 (blue circles) and 91 (red 
crosses). In addition, LDA predicts small hole pockets due to band 
93 (green rhombi). Each band is doubly degenerate in the spin-orbit 
mixed basis, hence only odd-numbered bands are shown. 

Fig. [3^)) contains two sheets (heavy bands 89 and 91) that 
are /-like in character, and a very small pocket (light band 
93, not shown) associated with boron p-electrons. The topol- 
ogy of the heavy sheets is non-trivial, featuring a combination 
of parallel quasi- ID sheets that are connected by a cylidri- 
cal quasi-2D tube. For comparison with future de Haas-van 
Alphen (dHvA) experiments we have determined the angular 
dependence of the extremal orbits 112211 . shown in Fig. Et) for 
the the lightest bands with with LDA-calculated m* < 10m e ). 
The heavier bands, with larger frequencies, are unlikely to be 
seen in dHvA measurements due to the quasi-particle damp- 
ing. 

Experimentally, the low temperature thermodynamics of (3- 
YbAlB4 is dominated by the spin entropy of the Yb ions||7|], 
indicating that this system is a Kondo lattice compound. We 
now determine the symmetry of the low lying magnetic dou- 
blets and construct the corresponding Kondo lattice model. 
In /3-YbAlB4 , the Yb ions are suspended between two seven- 
member rings of boron atoms (Fig.Q]) giving rise to an approx- 
imate 7 -fold symmetry. In this low symmetry environment, 
we expect the eight states of the J = 7/2 multiplet to split 
into four Kramers doublets. The approximate local symmetry 
group of the Yb ions is Cj x C^v . The only crystal field opera- 
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tors with both 2 tt/ 7 rotational and time reversal symmetry are 
polynomial functions P[J%] of J% 112311 . so to good approxi- 
mation, the magnetic doublets are eigenstates of definite mj. 
Empirically, the observation of a Schottky peak in the specific 
heat at Ta ~ 30-40 K, indicates that the lowest lying crystal 
field excitation has energy A w 2/c^Ta ~ 60-80 K. 

To determine the symmetry of the lowest lying Yb dou- 
blet, we appeal to geometric considerations. The lowest en- 
ergy magnetic doublet will correspond to the 4/ wavefunc- 
tion with maximum overlap with the boron rings. This rules 
out the states with \mj\ = 1/2 which are aligned along the 
c-axis, and the states with \mj\ = 7/2 which are planar, lying 
in the ab plane. The remaining contenders with mj = 3/2 
or 5/2 are "conical dumbbells" aligned along the c-axis. The 
seven-fold boron rings subtend an angle of 50° about the c- 
axis, closely corresponding to the aperture angle (52°) of the 
conical wavefunction of the | ± 5/2) multiplet. We expect this 
state to have the maximum hybridization with the boron rings. 
The most likely crystal field excitation state is the mixture of 
| ± 1/2) and | ±3/2) configuration, with the second-largest 
overlap. This simple model of the Yb ion provides a good 
fit to the experimentally measured magnetic susceptibility at 
high temperatures (Fig.©, provided one fits the Curie-Weiss 
temperature 0* w ~ —230 K due to antiferromagnetic RKKY 
interaction between Yb moments. The best fit is obtained with 
the excited state being mainly | ± 1/2), with a small admix- 
ture of the m j = ±3/2 as indicated in the inset of Fig. ffl with 
7 ~ 0.28. This is likely due to crystal fields of Al atoms that 
break the seven-fold symmetry of boron rings. 

Given this single-ion picture, the simplest effective model 
for the low-energy degrees of freedom is the one in which low 
lying | J = 7/2, m j = ±5/2) 4/ doublet interacts with the 
neighboring layers of two-dimensional boron p-electrons via 
a Kondo exchange interaction. Based on our calculations, we 
expect the spins on neighboring Yb ions to interact via an in- 
plane antiferromagnetic RKKY interaction . This leads us 
to propose a "layered Kondo lattice model" for /3-YbAlB4 , 
H = H ± H K + H H , where 

p=±, {m,n)e(ab) 

H K = JkSj I Yl C ]ap Cr a(3C j f3p+ Yl C ]ap a <xf3Cjf3p> 
\j,p=± j.P^p' 

Hh = ^ JijSi • Sj (2) 

(ij)e(ab) 

where the cj.(r mp ) operators create conduction electrons with 
spin a = ±1/2 in the two dimensional boron planes at site 
r m ± = r m ± c/2. The Sj are the Yb spin operators, interact- 
ing with boron layers via the Kondo coupling Jk- The {cj ap } 
represent boron Wannier states in the p th layer projected onto 
Yb site Rj in the \a = ±|) state, given by 
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FIG. 4: (a) Theoretical fit to experimental susceptibility data 0] 
Xc(T) (circles) and Xab(T) (rhombi) in the range 50-300 K us- 
ing the single-ion model shown in the inset and the fitted Curie- 
Weiss temperature 0cw- The best fit yields the lowest crystal field 
splitting A = 80 K, the coefficient of admixture 7 w 0.28 and 
#cw — —230 K, close to the experimental value 0^ ~ — 210 K. 



where summation is over atomic positions of the seven- 
member boron ring above or below the Yb site. Here the form- 
factors y h Ja (f) = c£ a Y lia - a (r), arise due to spin-orbit cou- 
pling [24], where Ca(a = (I = 3, a — a; a\ J = ~, a) are the 
Clebsch-Gordan coefficients for the J = 7/2, / = 3 state and 
Yi m (r) are spherical harmonics. Written explicitly 



= J_f \/6F 3 ,2(r) 
V7 V Y 3 ,- 3 (r) 
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The angular probability distribution function of the | ± f ) 

doublet is given by ^Z aa |3^a 2 (r) | 2 as shown in Fig.[T] 

The above model © provides a starting point for future 
theoretical studies of /3-YbAlB4 . Here we focus on the qual- 
itative physics of the model. At high temperatures, local mo- 
ments form on the Yb ions, which interact via RKKY inter- 
action, whose scale is however too small to allow magnetic 
ordering. The high temperature electronic behavior is then 
governed mainly by the "small" Fermi surface formed out 
of the quasi-two-dimensional boron bands, coupled weakly 
along the c-axis via the interplanar aluminum atoms. As the 
temperature is lowered below the Kondo temperature Tk ~ 
jj e -D/2j K (h ere ]j i s th e conduction electron band width), 
the quenched local moments will give rise to resonant bound- 
states that link neighboring boron layers, forming a three- 
dimensional Fermi surface of heavy quasiparticles. 

The observation of field-tuned quantum criticality in this 
material indicates that the RKKY interactions and the Kondo 
temperature are of comparable size, estimated [7] to be of 
order of 100 K, a value substantially larger than other quan- 
tum critical systems, such as YbRfi2Si2 B25I1 (Tk ~ 25 K) or 
CeCue-^Au^ [2] (Tk ~6K). The enhanced Kondo temper- 
ature suggests that this system balances a more mixed valent 
Yb ion with stronger RKKY interactions. Set against these 
observations, the tiny size of the superconducting transition 
temperature is quite mysterious, for a large magnetic interac- 
tion might be expected to induce stronger pairing |5|]. 

A striking feature of this compound is its highly 
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anisotropic, low- symmetry structure. The layered Kondo 
model Eq. © is akin to a quantum dot symmetrically cou- 
pled to two leads 12611 . where at high temperatures or in an 
applied field the /-electrons are localized and the system is ef- 
fectively Coulomb-blockaded. As the Kondo effect develops 
at low temperature, the interlayer channels opened up by the 
Yb Kondo resonances will lead to co-tunneling processes de- 
scribed by the second term in Hk, enhancing the c-axis con- 
ductivity. In /?-YbAlB4 , neighboring boron planes are weakly 
coupled by the s,p electrons of the Al and Yb atoms, leading 
to a finite c-axis conductivity even in the absence of the Kondo 
coupling. We expect however that the interlayer Kondo effect 
will resonantly "short-circuit" the boron planes at low temper- 
atures. In the light of the above discussion, the conventional 
scenario of the quantum critical point (QJ would suggest an 
immediate drop in c-axis conductivity as the system is tuned 
into the magnetic phase by applied pressure. 

The most puzzling feature of this material is the quantum 
critical behavior observed under ambient conditions. It is con- 
ceivable that the system is fine-tuned to a quantum critical 
point. In this case, pressure tuning will immediately drive (3- 
YbAlB4 into an antiferromagnetic phase, with a mass diver- 
gence following the form Eq. (OD), but where B C (P) is now 
finite. In this scenario, the fact that B c = at P = would 
be purely accidental. A more intriguing possibility, is that the 
field-tuned properties of /3-YbAlB4 are signature of a robust 
quantum critical phase. In this case, the critical field would 
be expected to remain zero over a finite range of pressure, 
while the exponent a may vary with pressure, according to 

m*(B)/m e oc B-^ p \ (5) 

The possibility of critical metallic phases in the Kondo lat- 
tice has recently been put forward by Anderson ll27ll . The 
existence of a paramagnet over a range of pressures between 
the fully developed Fermi liquid and antiferromagnet has also 
been proposed in the context of "Kondo breakdown" models 
of quantum criticality ll28[l29ll . however in that case the inter- 
mediate phase is a non-critical spin liquid phase. 

The existence of quantum critical point at zero field is ac- 
tually well known in the context of the two impurity Kondo 
model EH], where a competition between the Kondo ef- 
fect and valence bond formation between the local moments 
leads to non-Fermi liquid behavior. Could this phase be stabi- 
lized by the formation of a spin liquid in the Yb layers? 

This work was funded by NSF grant DMR 0605935. The 
authors are grateful to Zachary Fisk for bringing Ref.[3 to our 
attention. A.H.N, thanks P.M.C. Rourke and S. R. Julian for 
kindly providing the authors with computer code for Fermi 
surface analysis. Note added in proof: after this Letter was 



submitted, a report on de Haas-van Alphen measurements by 
O'Farrell et al. I32I1 has become available, in good agreement 
with our prediction of two heavy sheets of the Fermi surface. 
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